March 21, 1997 UMN-TH-1531/97 
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N flavor QED in two dimensions is reduced to a quantum mechanics problem with 
N degrees of freedom for which the potential is determined by the ground state of 
the problem itself. The chiral condensate is determined at all values of temperature, 
fermion masses, and the 8 parameter. In the single flavor case, the anomalous mass 
(m) dependence of the chiral condensate at 8 = n at low temperature is found. 
The critical value is given by m c ~ .437 • e/ifn. 



1 QED 2 and Quantum Mechanics 

Two-dimensional QED with massive fermions is not exactly solvable. The 
relevant parameter measuring the strength of the interaction is e/m where e 
and m are the gauge coupling and fermion mass, respectively. The massless 
limit corresponds to the strong coupling, whereas the large fermion mass limit 
corresponds to the weak coupling. In both limits the model is exactly solv- 
able. For the intermediate coupling e/m = 0(1), however, the model is highly 
interacting and some approximation methods need to be employed. 

The model has many similarities to the four-dimensional QCD in such 
respects as chiral dynamics and confinement. In this article we would like to 
focus on the chiral condensate at all values of the coupling, temperature, and 
6 parameter. 

The Lagrangian is 

N 

£ = 4F F F" + ^^{ 7 ' , (!a M - ei,) - m a }i . (1) 



=1 



We examine the model defined on a circle S 1 with a circumference L. Upon 
imposing periodic and anti-periodic boundary conditions on the bosonic and 
fermionic fields, respectively, the model is mathematically equivalent to a the- 
ory defined on a line (R 1 ) at finite temperature (T) by analytic continuation 
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to imaginary time t (= it) and the interchange of t and x. Various physical 
quantities at T 7^ on R 1 are obtained from the corresponding ones at T = 
on S 1 by substituting T -1 for L. 

In a series of papers it has been shown that the field theory system (|l|) is ef- 
fectively reduced to a quantum mechanical system of N degrees of freedom. 1-5 
The argument proceeds as follows. 

Fermions are bosonized on a circle. Each two-component fermion (ip a ) is 
expressed in terms of zero modes (qt'Pa) an d oscillatory modes (<p a (x),H a (x)). 
The boundary conditions enforce that p^s take integer eigenvalues. 

The relevant parts of the Hamiltonian are expressed in terms of {q a =Qa + 
1aiPa=\[pi +PZ])i (^a^a), and (&w,Pw) where 9 W is the Wilson line 
phase, the only physical degree of freedom associated with gauge fields on a 
circle. In the massless fermion theory the zero modes and oscillatory modes 
decouple. There appear one massive oscillatory mode, xi with fi 2 = Ne 2 /it, 
and N — 1 massless oscillatory modes, \a with fi^ = (a = 2 ~ N). 

Fermion masses introduce nontrivial interactions among these variables. 
All of the oscillatory modes become massive. Fermion masses change the vac- 
uum structure. The chiral condensate is induced and the boson mass spectrum 
is modified. Each affects the other, and must be determined self-consistently. 

There are N + 1 relevant zero modes, and q a . The 9 vacuum struc- 
ture, which follows from the gauge invariance, eliminates one degree. With an 
appropriate choice of the basis, we have N zero mode degrees which we denote 
by Pw and (p a (a = 1 ~ N — 1). The vacuum wave function must solve the 
Schrodinger equation 



H f{pw,<p) = ef(pw,<p) 

d 2 _ 47T 2 (iV 

^2 



d 2 4ir 2 (N~l) 
H = -l&2 — " \ + V N (p w ,<p) (2) 



where 



d 2 2 ^ d 2 



1 dip 2 ~ N-l ^ d^ a 
a=l r " a<b 



/TTfiLp w \ 2 4tt^ - / 2irp w \ 
V N = {^^) -■ w l^m a LB a coB[<p a -—) 

0=1 

N-l 

= cS ~ ^2 ip a . (3) 

a=l 

In the potential Vjy, B a depends on the boson mass spectrum fi a as well as 
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e, L, and m a . The detailed form has been given in refs. 1 and 3. The spectrum 
fi a is to be determined from the vacuum wave function solving (|^). We have 
a routine 

V(p w , <fi) -> f(pw, < jf) > A*q * V(p w ,ip) . (4) 

This is a rather non-trivial condition to be satisfied. The Hamiltonian is de- 
termined by its ground state wave function. 



2 Chiral condensates — T, m, and ^-dependence 

In the strong (e/m 3> 1) and weak (e/m <C 1) coupling limits Eq. (^) and the 
condition (^) can be solved analytically. For general values of the parameters, 
however, they must be solved numerically. We have numerically solved the 
Schrodinger equation on workstations. 

Let us focus on the N = 1 (single flavor) model. In this case there is no cp 
degree, the Schodinger equation (g) becoming 

d 2 1 

- ~rr + V (Pw) > f(pw) = e f{pw) 

V{pw) = (nuLpw) 2 — k cos(8 — 2-Kpw) 
k = 4:TTmLB(fj,iL) 

. . Z ( 7T f°° d,U \ 

B(z) = — exp 7 + - - 2 / 5 

where /i 2 = e 2 /7r. Note that B(0) = 1. /ii is the boson mass which is deter- 
mined by the vacuum wave function f(pw)'- 

, o 8ttitiB(uiL) . , , . 

Mi = A* + ( cos(27rp w - 0) )/ . (6) 

In the weak coupling limit e/m <C 1 at T = (L — oo), ( cos(27rpw — )/ = 1 
and the boson mass /zi = 2e 7 m. On the other hand, in the strong coupling limit 
e/m 3> 1 or in the massless fermion theory, ( cos(2npw — &)) f = e~^lv- L cos#. 
Eqs. (||) and (||) must be solved simultaneously. 
The chiral condensate is given by 

{^^) e = - 2B ^ lL) ( cos(27rp w -6)) f + —m. (7) 

The condensate is normalized such that it vanishes in the weak coupling limit 
at T = 0. (In the earlier references 1-5 the chiral condensate has been defined 
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without the last term in (Q). This reflects ambiguity in defining composite 
operators.) In the massless theory (m = 0), /ii = /i and 

(i>i/))0 = -yB(fj,L)e- v ^ L cos0 (8) 
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Figure 1: T dependence of the chiral condensate in the N = 1 model. Lines and points are 
for 6 = and tt, respectively. The mass m in the figure is measured in the unit of fi. 

In fig. 1 the condensate (tp if)) g j ' p<, is plotted as a function of the temperature 
T/fi = 1/fiL. There is a crossover transition around T/^i ~ 1. At high 
temperature the condensate becomes 9 independent, irrespective of m/fi. At 
low temperature, however, there appears a significant difference between 9 = 
and 9 = tt. 

At low T the condensate approaches — fj, cos 9/ 2tt in the massless theory. 
It, however, vanishes in the large mass limit as it approaches a free theory. 

At high T the condensate vanishes in the massless theory, whereas it ap- 
proaches a non-vanishing value in the massive theory. Thermal excitations 
yied finite condensates in the massive theory, as there is no chiral symmetry. 

The 9 dependence of the condensate originates from the potential. At low 
T//i = 1/fiL <gc 1 it is given by 

V = ir 2 ii 2 L 2 p 2 w - m/iiL 2 cos(6> - 2tt Pw ) . (9) 

At 9 = the potential is always minimized at pw = and everything smoothly 

changes. As the fermion mass becomes larger, the condensate increases. 
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Figure 2: m dependence of the chiral condensate in the N = 1 model. T/fi = .03 

On the other hand at 9 — -k there appear two degenetrate minima for 
large to. The problem becomes delicate for a moderate value of m/fi ~ 1. /ii is 
determined by (|^) which involves the vacuum wave function f(pw) determined 
by the potential (Q) itself. There appear two or three consistent solutions. 

For small m/n, the condensate decreases as T/fi increases. The behavior 
is opposite for large m/fi. The condensate increases as T/fi. For intermediate 
m/fi, the condensate initially increases as T/fi, but eventually starts to de- 
crease. We recognize profound structure in the behavior of the condensate for 
the intermediate coupling e/m =0(1) at 9 = it. 

3 Anomalous behavior at 9 = it 

To get more insight into the anomalous behavior of the condensate near 9 = ir, 
we have depicted the to dependence of the condensate at T//i = .03 in fig. 2. 

At 9 = the absolute value of the chiral condensate at T = decreases 
as m/fj, increases, and vanishes at m/fi — oo. This result has been, previously 
obtained by Tomachi and Fujita by the Bogoliubov transformation^ Our result 
agrees with theirs numerically. 

However, at 9 — n, there appears a discontinuity in the to/ fj, dependence. 
To magnify this point, more detailed study is presented in fig. 3, in which the 
condensate is plotted for 0.4 < m/fj, < .48 at T/fi = .03. The discontinuity 
occurs at m/fj, = .437. 

Why and how can such a discontinuity arise in a theory with a Hamiltonian 
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Figure 3: Discontinuity in the m dependence of the chiral condensate in the TV = 1 model 
at 8 = 7T. 



having the smooth dependence on m? To understand it, we have to go back 
to the equation (||). 

Suppose that fj,L = fi/T, m/fi, and 9 are given. The potential V(pw) in 
(||) is then fixed by the coefficient k of the cos(0 — 2irpw) term. With a given 
«in, the vacuum wave function f(pw) is determined, solving the Schrodinger 
equation. Now the boson mass fi\ (or [i\L numerically) is obtained by solving 
@ from f(pw)- The output K out = 4:TrmLB(fiiL) is determined with this new 

This process defines a mapping K out = g(Ki n ): 

Kin -> V(pw) -* /(Pw) -* Ml£ -> Kout • (10) 

We are looking for a solution K out = Ki n , namely a fixed point of g(n). 
Although the function g{n) depends on m//i very smoothly, there takes place 
bifurcation in the fixed point structure as m/fi varies. 

This problem was investigated in ref. 5. It has been found that for T/jj, < 
0.12 and 9 — 7r, there appears two attracting and one repelling fixed points 
for m c J fx < m/fi < m' c //i. Among these three fixed points, the biggest k 
corresponds to the lowest energy density and is chosen. Hence the lower critical 
mass m c / fx appears as the location of the discontinuity in the mj \x dependence 
of the various quantities. In fig. 4 we have displayed the m dependence of the 
k parameter at T / fj, — 0.03. 

In general m c depends on T. At T = it can be determined analytically 
to be m c //i = 0.435. At finite T we have found numerically m c //i=0.437 and 
0.454 at T/|U=0.03 and 0.07, repectively. 
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Figure 4: Discontinuity in the m dependence of the k parameter in the N = 1 model at 
6 = it. 

This gives us a puzzle. The Mermin- Wagner theorem ensures that there is 
no discontinuity in the T dependence of any physical quantities. It implies that 
if there is a discontinuity in the m dependence at m c , m c must be independent 
of T. In our approximation we have found that m c is almost universal, but has 
small T dependence. Indeed, we see a jump in the T dependence of the chiral 
condensate for m/fi = .437 in fig. 1, though it may be a smooth transition 
when the approximation is improved. 

There are two possible scenarios. The first possibility is that m c is uni- 
versal and T independent in the full theory. The second possibility is that the 
discontinuity in m disappears and is replaced by a rapid crossover in the full 
theory. It is a challenge to know which picture is real. 

4 Generalization and summary 

The generalization to the N flavor case is straightforward. Extensive analysis 
for small fermion masses has been given in refs. 1-4. For instance, the chiral 
condensate at low T for degenerate fermion masses (m a = m /i) is given by 

1 .— 1 / V \ Tv+T /m \ Tv+T 1 

— {iptp/e = 2e 7 cos — — forT«m»+ifi»+i (11) 

fjj 4ir \ AT / v n J 

where 8 is defined in the interval — it < 9 < +ir by 6 = 9 — 2tt[(9 + ir)/2ir}. 
Recently SmilgaQ has obtained an exact result for N = 2 at 9 = 0, T = 0, and 
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m <C /i, which agrees with our result within 5 %. At T = the cusp singularity 
appears at 9 = tt. The A = 3 model with fermion masses mi < mi -C TO3 
mimics the four-dimensional QCD. 

In this article we have presented an alternative method to explore QED 
in two dimensions. QED is reduced to a quantum mechanical system of fi- 
nite degrees of freedom whose Hamiltonian needs to be determined by its 
ground state wave function itself. The associated Schrodinger equation has 
been solved numerically on workstations. With this algorithm one can deter- 
mine various phyisical quantities at any temperature and 6 and with arbitrary 
fermion masses. Our method supplements other numerical methods such as 
the lattice gauge theory and light-front quantization method. 8-9 More results 
by these methods are wanted for comparison. 
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